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Abstract. In terms of the crystal base of a quantum affine algebra Uq{g), 
we study a soliton cellular automaton (SCA) associated with the exceptional 
affine Lie algebra g = ' . The solitons therein are labeled by the crystals 
of quantum affine algebra Uq{A^j^^). The scatteing rule is identified with the 
combinatorial R matrix for C/q(A^'^')-crystals. Remarkably, the phase shifts in 
our SCA are given by 3-times of those in the well-known box-ball system. 



1. Introduction 

1.1. Box ball systems in terms of crystals. The box-ball system [19, 21] is 
a discrete dynamical system in which finitely many balls move along the one- 
dimensional array of boxes under a certain rule. Then each isolated array com- 
posed of balls behaves like a soliton: (i) The longer isolated soliton moves faster, 
(ii) The number of solitons does not change under the time evolution. (Hi) If the 
solitons are apart enough from each other at an initial state, then their lengths do 
not change after their collisions. Its integrability has been proved by making a 
connection to the difference analog of the Lotka-Volterra equation through the lim- 
iting procedure called ultra-discretization [22]. In the last part of 90's, the box-ball 
system was re-formulated [8] in terms of crystal base theory [9, 14, 15] associated 
with the symmetric tensor representations of the quantum affine algebra Uq{g) for 

non-exceptional affine Lie algebra 0„ = Bn \ Cn \ Dn \ ^2ri-i' ^^^^ ^i+i- 
Let us recall the box-ball system in terms of crystals. For simplicity we restrict the 
discussion in 0„ = An '' [2, 3]. By a similar discussion, there are some expansions 
[7, 24]. 

Recall necessary results on box ball system in terms of crystals. For Uq{An^)- 
crystal Bi {I G Z>i), see section [2^ Fix a sufficiently large positive integer L. Put 
ui — {1,0, . . . ,0) e Bi. Define a set of paths Tl by 

^ { P = bi(Eib2'S> ■ ■■ <»bL e Bf^; = ui if j > 1 } . 

By the iterating crystal isomorphism Bi >Si Bi ~ i?i (g) S;, we have a map from 
Bi ® ?L to Tl ^ Bi by 

Ul ® &i®&2 ® • • • ® &L 

^bi®u'i®b2®---®bL^ ^bi®b2® ■ ■ ■ ®bL®ui. 

Thus we define a map Ti : CP^ — > CP^ by 



Ti{bi ®b2®---®bL) = bi®b2®---®bL. 
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Application of Ti induces a transition of state. Thus map Ti can be regarded as a 
dynamical system, in which T; plays the role of time evolution. We call the system 

{T/(p) I t > 0, ; G 



1, p e ^l} the An'' -automaton. Thanks to the Yang-Baxter 



equation, {Tj}j>i becomes a commuting family. 

A m-soliton state of length {li,l2, . . . ,lm) is given by 



...[h] [h] ••• [U-- 

Here . . .[I]. . . denotes a local configuration such as 
(1.1) • • • (8) (1) O (1) ® {n + I)"' " X ■ ■ ■ : (.S)'"'-''^ 



Xl+X2-\ \-X„ 



sandwiched by sufficiently many (l)'s. Here (■) signifies a frame of Yoimg tableau. 

The state Tk{p) is obtained by the rightward shift by min(fc, I) lattice steps. The 
following figure explains the action T; {I = 3) on a one-soliton state (It transfers 
the upper line into the lower line): 
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The array "332" behaves as a soliton just like in the classical soliton theory, where 
the frames are omitted. 

Let Bi {I e Z>i) be a i7^(A^^2i)-crystal. Define a map h : BiU {0} Bf U {0} 

by 
(1.2) 



(n -f 1)®="" (g) • • • (g) (3)®="2 ^ (2)®^i if 6 7^ 0, 
if 6 = 0. 



Then it is easy to sec that for any ie{l,2, ...,n — 1} and 6 G i?; U {0} 

(1.3) Liie,{b)) ^ e,+i{ii{b)), niUb)) ^ h+i{ii{h)) 

Thus local states ...[/]... are labeled by crystals of the smaller algebra Uq{A^^}_-^). 
We also denote m-soliton state with length (Zi,/2, ■■■■,1m) by the tensor product of 
f^9(^n-i)-crystals. 

(1.4) z'^^hi (8> z'^^b^ (8> • • • O ^^""6™ G Afr(S;J Aff(BiJ • • • hS{BiJ 

Here 7 denotes a phase of the local configuration ...[I]... and defined by 7 = 
min(r, Z)t -|- (the position of "...[Z]..." from the right end) under T^. It is invariant 
unless ...[Z]... interacts with other solitons. 

For the m-soliton state with length (^1,^2, ■ • ■ , ^m)) we assume li > I2 > ■ ■ • > Im^ 
Since the longer solitons move faster under the time evolution (r > Zi), we can 
expect that the state turns out to be 



(1.5) 



® z''% (g) G AS{BiJ g) • • • ® Aff(B(J (g) AS{Bi,). 



after sufficiently many time evolutions. We represent such a scattering process as 
follows. 



(1.6) 
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Theorem 1.1. ([2, 3]) For m > 2 the map Sm is factorized into 82- The map 
S2 is described by the combinatorial R for Uq{A^^}_^)-crystals. Especially, the phase 
shift is given by 

(1.7) 7^-72=71-7^ = 2^2 + i^(6l®62)• 
1.2. Our case. In [16] a family of J7^(£'4^-')-crystals {i?;}/>i corresponding to the 
vertex 1 in the Dynkin diagram was studied. Kashiwara, Misra, Okado and the 
author gave an exphcit description of the crystal structure by the coordinate repre- 
sentation (see. section [2. 3p in the spirit of [11]. Thus it is expected that a soliton 
cellular automaton associated with J7^(D4'^'')-crystals can be constructed by the 
same method as [3, 8]. 

First, we have to make a table of a combinatorial R matrix for [/^(D4'^'')-crystals. 

R : AS{Bi) ^ Aff(Si) — > Aff(Bi) ® Aff(S/) 

Now we introduce an algorithm for the combinatorial R (section [3|) by the Lecou- 
vey's column insertion algorithm for type G2 [17]. We describe it only for those 
tableau with shape (n) or (n + 1, 1) for n G Z>o, which suffices for our aim. The 
proof goes as follows. As C/g(G2)-crystals, Bi ®Bi decomposes into connected com- 
ponents that are isomorphic to the crystals of irreducible J7q(G2)-modulcs. Within 
each component, the general elements are obtained by applying /j's (i = 1, 2) to the 
Uq{G2) highest weight elements. Thus we first prove the theorem directly for the 
highest elements (See Appendix A). For the general elements, the theorem follows 
from the fact due to [17]. (See Proposition 13.91 and Lemma [3. 101 ) 

Next, using the algorithm of combinatorial R given in the previous section, we 
construct a soliton cellular automaton associated with Ug{D\ )-crystals. We call 
it the -04^' -automaton. The method is same as the An"^ case. We summarize the 
result below: (We use the same notation as the aIP case.) 

• A one-soliton state in D\ -automaton is described as 

(1.8) ■ • ■ ® (1) ® (1) ® (3)^"^ ® (2)»"i (E)(1) ® • • • ^ z^{xi,X2) e Aff(^0. 

Xi+X2 = l 

Here Bi denotes the U^{A^^'')-CTystal. 

• The two-body scattering of solitons is described by the combinatorial R for 
[/^(A^^^)-crystals. Especially, the phase shift is given by 

(1.9) 72 - 72 = 71 - 7i = 2/2 + 3 X H{bi ® 62)- 

• Scattering of solitons is factorized into two-body scattering. 

Let us compare the Ag"'^'- automaton with _D4^''-automaton. Let {Bi}i>i be the 
family of [/^(yl4^'')-crystals. Under the time evolution T2, we give an example of 
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two-soliton scattering described by AS{B2) ® Aff(Bi) Aff(Bi) Aff(S2). 



= 
= 1 
= 2 
= 3 
= 4 
= 5 
= 6 
= 7 
= 8 
= 9 



t 
t 
t 
t 
t 
t 
t 
t 
t 
t 

t = 10 
t = 11 



1)4'^ -automaton 

22111311111111111111111111 
11221131111111111111111111 
11112213111111111111111111 
11111122311111111111111111 
111111112Q1111111111111111 
1111111111|111111111111111 
11111111111^21111111111111 
11111111111110211111111111 
11111111111111232111111111 
11111111111111121321111111 
11111111111111112113211111 
11111111111111111211132111 



^2^^ -automaton 



221113111111111111111111111 
112211311111111111111111111 
111122131111111111111111111 
111111223111111111111111111 
111111112321111111111111111 
111111111213211111111111111 
111111111121132111111111111 
111111111112111321111111111 
111111111111211113211111111 
111111111111121111132111111 
111111111111112111111321111 
111111111111111211111113211 



Here the markers specify the positions of initial sohtons under the interaction-free 
propagation. We can read off the phase shift from the deviation of the outgoing 
solitons from the corresponding markers. In terms of the soHton labels, the above 
scattering is described by the combinatorial R for U q{A^p )-CT:ysid\: 



Aff(B2) Aff(Si) 
j^-2(2,0)®z^-^(0,l) 



Aff(i?i) ® Aff(i32 
(l,0)®z 



(L-6)+5 



(L-2)-5 



(1,1)- 



Here the phase shift 5 is given by 

5=2x l + i/((2,0)® (0, 1)) X 



3 if = 7^ 
1 if a = a: 



(3) 
4 ' 
(1) 



1 if I 



D 



(3). 



1 if 



A. 



(1) 



In order to prove the theorem, we introduce an operator T<^ constructed from the 
C/^(D4'^^)-crystal which corresponds to the vertex 2 in the Dynkin diagram. The 



idea is same as [7]. is to put the highest element y^j G on the left and to 

carry it through to the right via the isomorphism B\^ ® Bf^ ~ Bf" ® B\^. It acts 
on the space of states and has the following properties: (See Example 14. 161 ) 

(1.10) Ti, commutes with T^. 

Z^{xi,X2) if X2 = 0, 



(1.11) 



T[,(z^(xi,X2)) 



z^-^{xi + l,X2-l) if 2:2 7^0. 



Here z^{xi,X2) denotes a one-soliton state labeled by the [/g(A^^^)-crystals. Since 62 
commutes with T^, we can reduce the proof to the C/q(Ai)-highest weight elements 

z^n',0) z^^(2/i,y2) e Afr(Bi) ® Afr(^fc) for / > k. 

From the second property of T\^, the proof is reduced to the induction hypothesis 
on the number of applications of . 
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1.3. Plan of the paper. This paper is organized as follows. In section 2 we recall 
the facts on the crystal and combinatorial R. We also recall a family of Ug{g)- 

crystals {i?(};>i for g = An'' or and C/^(_D4^-')-crystal associated with the 
vertex 1 and 2 in the Dynkin diagram respectively. In section 3, we construct an 
algorithm for the combinatorial R matrix in terms of Lecouvey's column insertion 
algorithm [17] for the tableaux of type G2- In section 4, using the algorithm in- 
troduced in the previous section we construct a box ball system associated with 
?7^(£)4"'^)-crystals and describe the scattering rule of solitons appearing in it. 
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2. Preliminaries 

2.1. Crystals. A crystal B is a set B with the map gj, /j : B U {0} — > B U {0}, 

satisfying the following properties: 

(2.1) g,0 = 7,0 = 0. 

(2.2) for any b and i there exist n > such that e"6 = /"6 = 0. 

(2.3) for b,b' eB and i, fib = b' b = e^b' . 

For 6 e B, set 

Ei{b) = max{n > | e^ft 7^ 0}, <^i(6) = max{n > | /?*6 ^ 0}. 

The crystal B is identified with a colored oriented graph, called the crystal graph, 
if one draws an arrow as 

(2.4) h^b' ^ fib = b'. 

For two crystals B and B' , the tensor product B (g) B' of is defined by as a set 
B ^ B' = {b^b' \ b G B, b' G B'} with the actions of e^, fi as follows: 

ei(6(g)6') 
Mb^b') 

Here 6 and 6' are understood to be 0. 

Let Q be an affine Lie algebra and let J7^(0) be the corresponding quantum affine 
algebra with out degree operator [14]. For J7g(g)-crystals B and B', assume that B 
(resp. B') is a finite set and B'S)B' is connected. Let Aff(i?) = {z% \ b e B, deZ} 
be the affinization of a t/g(fl)- crystal B. Then Aff(B) is an infinite set, and is 
endowed with iiJi by ei{z%) = z'^+^''' {iib) , fi{z% = z'^-^'°{fib). 



e^b®b' li^^{b)>e^{b'), 

b(E)e^b' if(pi{b) < e^ib'). 

hb®b' if^,(6) >e,(fe'), 

b®~f,b' ii^,{b)<e,{b'). 
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The combinatorial i?-niatrix (combinatorial R for short) is a map with the fol- 
lowing form: 

R : Afr(B) (g) Aff (B') — > A&{B')^m{B) 

Here a map (. : 6 (g) 6' i— > 6' 6 is a crystal isomorphism, manely it commutes with 

Cj, fi. H{b 6') is a Z- valued function on B ® B' satisfying the following property: 
for any b®b' & B®B' and i € I such that ei{b (g) b') = 0, 

r JJ(6 ® 6') + 1 if i = 0, (po(fc) > eo(&'), Mb') > £o{b) , 
H{ei{b ® 6')) = < 6') - 1 if i = 0, ipo{b) < eo{b'), Mb') < eo(b), 
[ H{b (g) &') otherwise. 

H{b (g) 6') is called the energy function, and is determined up to a global additive 
constant. The existence of the crystal isomorphism i : B(g)B' — * B'(gB and energy 
function H{b (g) b') are guaranteed by the existence of the R matrix in [14]. The 
definition of energy function assures that the combinatorial R commutes with Cj, 
Denote the crystal isomorphism 

B (g) B' B' (g) B, b®b' ®b 

by the following symbol. 

b' 



b' 

2.2. J7^(A^^-')-crystal Bi. Let {i?;};>i be a family of crystals of the Z-fold symmet- 
ric fusion of the vector representation of U'q{An^). It is as a set 

-Bi = = {xi,X2, ■ . ■,Xn+l) I Xi e Z>o, = /| . 

The actions e,, /j (i = 0, 1, . . . , n) are given by 
eob ={xi - 1, . . .,Xn+i + 1), 

Bib ={XI, ...,Xi + l,Xi+i - 1, . . .,Xn+l) if i 7^ 0, 
fob ={xi + 1, ....... , Xn+l - 1), 

fib ={xi, ...,Xi- l,Xi+i + 1, . . .,Xn+i) if i 7^ 0. 

Let Aff(B;) be the afBnization of B;. A piecewise linear formula [3] for the 
combinatorial R matrix for a crystal Bi Bi> 

R : AS{Bi) ^ AS{Bi,) AS{Bi,) ^ AS{Bi) 

is given as follows: 



SCATTERING RULE IN SOLITON CELLULAR AUTOMATON ASSOCIATED WITH CRYSTAL BASE OF (/.((D^^'l 

for 6 (g) 6' = {xi,X2, ■ ■ . ,a;„+i) ® (yi,2/2, ■ • ■ ,2/n+i) & Bi® Bi,, 

{k-l n+l 1 

^ — ' ^ — ' I 
J=l j=k+l ) 

y ®h= {y[,y2, ■ ■ ■,yn+i) ® {x[,x'^, . . . 

x[^x, + Q,{b, b') - Q,.i{b, b'), y'i = yi + Qi-i{b, b') - Qi{b, b'), 
H{bi®b2) = -Qn+i{b, b'), 

where Xn+i+j (resp. yn+i+j) denotes Xj (resp. yj) for 1 < j < n + 1. 

Another algorithm is known in [18] , which is described by the tableau represen- 
tation of Bi . Then, the element b £ Bi is described by the Young tableau of shape 
(Z) with entries 1, 2, . . . , n + 1. 

(2.5) {xi,X2,...,Xn+i) € Bi (11...1122...22 ■■■ n + 1 n + 1 ) 

Xl X2 Xn+1 

Here the frame of tableau is replaced by ( • • • ). 

Example 2.1. The n = 3, I = A, I' = 2 case. 

• Coordinate representation. 

z'^{2, 1, 1, 0) O ^'''(0, 1, 1, 0) ^ ^<*'+°(l, 1, 0, 0) O ^''-°(1, 1, 2, 0), 
z'^{2, 1, 1, 0) (g) z'^'il, 1, 0, 0) ^ z'^'+'^il, 0, 1, 0) O z'^-'^{2, 2, 0, 0), 
z'^iO, 1, 1, 2) O z'^'il, 1, 0, 0) ^ z'^'+^{0, 0, 0, 2) O z'^-^{l, 2, 1, 0). 

• Tableau representation. 

z'^(1123) ® z'^'{2:i) ^ z'^'+i(12) O z'^-\U3'S), 

z\n2S) (g) z'''(i2) ^ z'^'+^n) (g) z'^-\n22), 

^'^(2344) ® ^'^'(12) z'^'+^iU) g) z'^-2(-i223). 
Now we can state 

Proposition 2.2. (Yang-Baxter equation). The following relation holds on Aff{Bi)® 
Ajf{Bv)(^Ajf{Br) : 

(2.6) {R g) 1)(1 (g> R){R g) 1) = (1 g) R){R (g> 1)(1 g) R). 

(3) 

2.3. U'q{D\ '^)-crystal Bi. A family of C/^(0)-crystals {-B(}z>i corresponding to the 

vertex 1 in the Dynkin diagram in the exceptional affine Lie algebra q — D^p was 
studied in [16] . In this paragraph let us recall the coordinate representation of the 
crystal B;. 

As a set crystal Bi is given by 

I 

(2.7) Si=0B(iAi) as C/g(G2)-crystals. 

3=0 

B{j Ai) is isomorphic to the crystal base for finite-dimensional irreducible Uq{G2)- 
module with highest weight jK\ in [12]. Any element of B{j Ai) is represented as a 
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one-row semistandard tableau with entries 1, 2, 3, 0, 3, 2, 1 as follows 
(2.8) 





2. ..2 


3-.-3 





3--.3 


2. ..2 






Wl 


W2 


W3 


Wo 


W3 


W2 





(2.9) 



Wo G {0, 1}, ^(Wi + Wi) +wo= j. 



Now we set 

(2.10) (a;i,a;2,a;3,i3,52,^i) = (wi,W2,2w3 + wo,2w3 + wq, w>2, w>i)- 

Define a coordinate representation of B(jAi) by 

{a;i,5i G Z>o, 
6 = (a;i,a;2,a;3,S3,52,Si) | X3 = X3 (mod 2), 
a;i + a;2 + (x3 + ^3)72 + X2 + xi = j 

Then the erystal structure is given as follows: 



eob ■■ 



fob- 



Xl 



1,...) 



if Zi + {Z2 + {3Z4 + {Z3)-)-)- < 0, 



, ,a;3 - l,a;3 - 1, . . . + 1) if 2:2 + {3z4, + {zi + Z3)-)- <0 < zi, 



.,X3-2,...,X2 + l,...) 
.,X2-l,...,X3 + 2,...) 
Xi - 1, . . . ,2:3 + 1,X3 + 1, . . .) 

■,Xi + 1) 

Xi + 1,...) 

. . . ,a;3 + 1,^3 + 1, . . . ,Si - 1) 
...,a;3 + 2, ...,^2 - 1,...) 

. . . ,.X2 + 1, • • • ,-^3 - 2, . . .) 
Xi + 1, . . .,X3 - 1,53 - 1, . . .) 

. ,Xl - 1) 



if 82:4 + {Z3 + {Zl)-)- <0<Z2 + {zi) + , 

if 3z4 + (z2 + (^1)+)+ > > Z3 + {zi)+, 

if Z3 + (3Z4 + {Zi + Z2)+)+ >0>Zi, 
if 2i + (z3 + (3z4+(z2) + ) + )+ >0. 

if Zl + {Z2 + (3^4 + (23)-)-)- < 0, 

if Z2 + {3Z4 + {Zl + Z3)-)- <0<Zi, 
if 3Z4 + {Z3 + {Zl)-)- <0<Z2 + {zi) + , 
if 3Z4 + {Z2 + {zl) + )+ >0> Z3 + {Zl)+, 

if Z3 + {3z4 + {zi + 22)+)+ > 0> Zl, 
if Zl + {Z3 + (3Z4 + (z2) + ) + )+ > 0. 



£0(6) =/ - s{b) + (zi + {z2 + 3z4 + {z3 + (zi) + )+)+)+ - (221 +Z2 + Z3 + 3Z4), 
(foib) =1 - S{b) + {Zi + {Z2 + 3Z4 + {Z3 + (21) + ) + ) + ) + . 



iib ■■ 



fib 



e2& : 



f2b 



(. . . ,X2 + l,Xi - 1) 
{...,X3 + 1,X3 - 1, . 
{xi + 1,2:2 - 1, . . .) 

(Xi - 1,X2 + 1, . . .) 



. ,X3 - 1,X3 + 1, . 
. ,X2 - l,Xi + 1) 

.,X3 + 2,X2-1,. 

. ,a;2 + l,a;3 - 2, . 



■,X2 
■,X3 



1, X3 + 2, 

2, X2 + 1, 



if X2 - X3 > (X2 - 2:3) + , 
.) if X2 - X3 < < X3 - X2, 
if (X2 - X3)+ < X2 - X3. 

if (X2 - X3)+ < X2 - X3, 
.) if X2 - X3 < < X3 - X2, 
if X2 - X3 > (X2 - X3) + . 

.) if X3 > X3, 
.) if X3 < X3. 

.) if X3 < X3, 
.) if X3 > X3. 
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ei{b) =xi + {x3 - X2 + {x2 - X3) + ) + , e2{b) ^ X2 + -{xs - 0:3)+, 

(pi{b) =Xi + (X3 - X2 + {X2 - ^3) + ) + , 'f2{h) ^ X2 + ^(X3 - Xz)+. 

Here (a;)+ = max(a;, 0), = min(a;,0), 

Zi=Xi-Xi, Z2^X2- X3, 23^X3- X2, Z4 = (^3 - X3)/2, 

X3 + X3 _ _ 

S(0) — Xi + X2 -\ \- X2 + Xi. 



For instance, as an element of i?'^^(7Ai) 



2 2 3 



(1,2,3,1,0,2). 



Let us denote the elements of Bi by 

□ = (1,0,0,0,0,0), 0= (0,1,0,0,0,0), 
[0]= (0,0,1,1,0,0), g= (0,0,0,2,0,0), 

|T|= (0,0,0,0,0,1), 0= (0,0,0,0,0,0). 

then, the crystal graph of Bi is given as follows: 



(0,0,2,0,0,0), 
(0,0,0,0,1,0), 




The arrows without number are 0-arrows. 

2.4. ;7^(Df ^)-crystal B^. It is known [6, 15, 16] that ;7^(Df ^)-module corre- 
sponding to the vertex 2 in the Dynkin diagram has a crystal pseudobase and is 
decomposed into ~ V{A2)®V{Ai)®'^ ®V{0) as highest-weight J7,(G2)-modules, 
where Ai = Ai — 2Ao and A2 = A2 — 3Ao. The crystal B^^ of V'^ is given by as a 
set B\f 

B^ =(A2) © (Ai)®2 e 5(0), 
BiO) ={0}, 

B(Ai)®2^ □ {(a), I a = 1,2,3,0,3,2,1}, 

3=1,2 



B{A2) - 




I The word a(3 is the element of -6(12). 



For a set of words i?(12), see section [STTl The labeling of the elements of B{A2) is 
taken from [12]. Thus the actions e^, fi {i — 1,2) on i?(A2) agree with those in [12]. 
On 5(0) and B(Ai)®2 they are defined by 

e,0 = 7,0 = (z = l,2), 

(1), ^ (2), ^ (3), i (0), ^ (3), ^ (2), i (1), U = 1, 2). 
Here h ^ b' means fib — b' . 
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The 0-action on is given as follows. 

(5)i(3,,ip).iQ. (5)i(2),i(3,.iQ, 

i) ^ (0), i (1)2, (0). i. @ , (l),i.*i.(l)„ 

°)^(3)., (5).i@. (°)i(2)., 

o)' (s)' (s)' (I 

2.5. Yang-BzLxter equation. For Ulj{Di )-crystals B we define the affinization 
of them by 

(2.11) AS{B) = {z%\d€Z,b€B}. 

Here z denotes spectral parameter. 

Define a combinatorial R matrix for a f7q(i?4^^)-crystals B (g) B' by 

R : AS{B) O Aff(B') — > AS{B') O Aff(B) 

z% ® z'^' b' ^ ^d'+Hib^b' ) I ^ ^d-H{b^b') y 

Thanks to the ordinary Yang-Baxter equation, we obtain the following result. 

Proposition 2.3. (Yang-Baxter equation). The following equation holds on Aff{B)® 
Aff{B')^Aff{B"): 

(2.12) {R 1)(1 (g) R){R (g 1) = (1 (g) R){R (g 1)(1 (g R). 
for B®B'®B" = Bi® Bi, ® Bi„ or B^^f^Bi® Bi, . 

Hereafter we normalize the value of energy function on Bi (g) Bii as follows. 

(2.13) H{{1, 0, 0, 0, 0, 0) ® (/', 0, 0, 0, 0, 0)) = 

3. Combinatorial R matrix 

3.1. Lecouvey's column insertion algorithm for type G2. We recall the col- 
umn insertion for the G2 case [17]. We define the set of words B{1), -B(IO), -8(12), 
i3(ll), _B(112) and B(121) as follows. Remark that a/37 denotes a word, and a /? 7 
an one-row tableau. 



B{1) = {1, 2, 3, 0, 3, 2, 1}, B(10) = {10, 13, 12, 22, 21, 31, 01}, 

5(12) = {12, 13, 23, 20, 23, 30, 33, 32, 00, 03, 02, 32, 31, 21}, 

/ 11, 21, 22, 31, 
^ > \ 21, 22, 23, 20, 



32, 33, 01, 02, 03, 31, 32, 33, 30, 33, 
23, 22, 11, 12, 13, 10, 13, 12, 11 
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' 112, 113, 212, 213, 223, 220, 223, 312, 313, 323, 320, ' 

323, 330, 333, 332, 012, 013, 023, 020, 023, 030, 033, 

032, 312, 313, 323, 320, 323, 330, 333, 332, 300, 303, 

302, 332, 331, 212, 213, 223, 220, 223, 230, 233, 232, ' 

200, 203, 202, 232, 231, 221, 112, ll3, 123, 120, 123, 

130, 133, 132, 100, 103, 102, 132, 131, 121 



B(121) 



121, 131, 122, 231, 201, 231, 232, 132, 133, 301, 331, 

332, 321, 322, 323, 232, 233, 001, 031, 032, 021, 022, 

023, 202, 203, 233, 230, 233, 321, 322, 323, 311, 312, 

313, 310, 313, 302, 303, 333, 330, 333, 320, 323, 322, 

211, 212, 213, 210, 213, 212, 002, 003, 033, 030, 033, 

020, 023, 022, 320, 323, 322, 312, 311, 211 



We also define the bijection ^ : B{10) B{1) and rj : B(121) B(112) by 
from the j-th component in B{10) (resp. B(121)) to that in B{1) (resp. 5(112)) 
respectively, e.g., ^(22) = 0, ^~H^)_ = 12, r/(021) = 030, Ty-^lOO) = 320. 

Given a letter a G {1,2,3,0,3,2,1,0} and at most two row tableau T, which 
contains 0, we define a tableau denoted by (a ^ T), and call such an algorithm a 
column insertion of a letter a into a tableau T. Let us begin with such T's that have 
at most one column. The procedure of the column insertion (a —>■ T) is summarized 
as follows. 

case (a ^ 0) = a 

a 



case 1 (/3 ^ a) = ^ if a/? € B(12) 
case 2 (/? ^ a) = (/3 a ^) ii a f3 e B {11) 



case 3 



case 5 (/?—>«) = < 



a 

7 ^ /? ; " V /3' 7' 
case 4 (/3 ^ a) = ( 7 ^ ); 7 = ^(a/3) if a/3 e B{10) 

P if a = (j), 
a if /? = 0, 
(0 if(a,/3) = (.^,0)or(l,l). 



; j'a'l3' = v{ctl3^) if a/37 G S(121) 



The above cases do not cover all the tableaux with two rows, but we only deal 
with these situations in this paper. When T is a tableau with shape (m, n) for 
n S {0,1}, we repeat tlie above procedure: we insert a box into the leftmost 
column of T according to the above formula. If it is not case or case 1 , replace 
the column by the right hand side of the formula. Otherwise, replace the column 
by the right hand side of the formula without the right box. We say that this right 
box is bumped. We insert it into the second column of T from the left and repeat 
the procedure above until we come to the case or case 1. 

Example 3.1. For any n G Z>o, let 61 18) &2 = {cti cti . . . a„) ® (/?) be the tableau 
representation of the elements in Bi ® By. Let h\ * hi be the tableau obtained by 
column insertion (62 — » Then we have 
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(3.1) 61 = ((/>) case: &i * &2 = if /? = </>; = /3 if /? 7^ 0- 

(3.2) 61 7^ ((^), 62 = (?!>) case: 61 * 62 = ai a2 • • • Q!„. 

(3.3) ai/3 G B(ll) case: 61 * 62 = /? ai a2 ••• 

(3.4) (ai,/J) = (1,1) case: 61 * 62 = Q!2 as ••• Q!„. 

(3.5) ai/3 e B{12) case: 61 * 62 = ^1 "2 as ... a„ 

(3.6) Q!i/? e -B(IO), Q!2C(ai/3) £ -B(ll) case: 61 * &2 = 7 0^2 • • . an- 

(3.7) aij3 G B(10), a2^(ai/3) e 5(12) case: 61 * 62 = ^2 ag a4 . . . a„ ^ 

where 7 = ^(ai/3) if ai/3 G B(10). 

We define the reverse bumping algorithm by RHS = LHS for the cases 0-3. 

Example 3.2. Fix n G Z>i. The procedure of bumping out all letters from the 
tableau 

Tn = K - (r'n-i - (• • • ir[ - {Qo - (po - 0))) •••))) 

is given by 

^ _ PnTn r„_i • • • r2 (ri ^ 0) 
"~ qn 

= Pnrnrr.-i---{ri+i^n)...n forl<i<n-l 

Qn ~ ~ 

Pn \ ( Pn-1 \ 

=(r; ^ T„_i) = (r; ^ (r;_i r„_2)) = • • • 
=(<->(<-i-(---W-ro)...))) 

=(r; ^ (r;_i ^ (• • • (r; ^ (go - (Po - 0))) •••))) 
where pi-\qi-ir[ = r]~^{riPiqi) for n > i > 1. 

3.2. Algorithm. Let {Bi}>i be the family of [/^(D4^'')-crystals. In this paragraph 
we will give the explicit description of the combinatorial R on Bi ® Bi. Fix the 
positive integer I. For a given element 61 (g) 62 in -B; (g) Bi, we apply the following 
procedure. 

Rule 3.3. Step 1. Denote 61 <8) &2 G (g) Bi by tableau representation 

(3.8) 61 (g) 62 = (ai Q!2 . . . a„) (g) (/?) G Si O Bi (1 < n < /)• 

Step 2. Compute the column insertion (62 -^bi), and obtain the tableau hi * 62- 
Let {p + q, q) be the shape of 61 * 62 . 

Step 3. Define a map if : Bi — *• Z by 



(3.9) iJ(6i (g) 62) = 

w/iere (61 * 62)1 = P + 9- 



-2 i/ai/3 e B(10), a27 G S(ll), 
max (—2, (61 * 62)1 — Z — 1) otherwise. 



SCATTERING RULE IN SOLITON CELLULAR AUTOMATON ASSOCIATED WITH CRYSTAL BASE OF U,(d'^% 



Step 4- Bump out all letters from 61*62 by the reverse bumping algorithm. Denote 
the procedure as follows: 

(3.10) 61 * 62 = (il ^ (t2 ^ (. . . {tp+2q ^ 0) 

Set (7r) = {tit2 ... tr); <r <p + 2q. 
Step 5. Define b'^ ® b[ € Bi ® Bi by 

{tp+2q) ® (Tp+2q_l) 
{tp+2q) ® {7p+2q-ltp^2q) 
(1) ® (Tp+2,1) 
(Tp+2g) <E) ((/)) 
(0) ® (Tp+29) 



(3.11) 



62 ® 6i = < 



•))); to = 



for case I, 
for case 11, 
for case III, 
for case IV, 
for case V, 
otherwise. 



whe 



III. 
V. 
I. 
II. 
I. 
I. 
III. 
IV. 

I. 

IV. 
V. 

I. 

V. 

III. 



aiP e B{11), n<l~l, 
aiP G 5(11), n = l-l, 
ai/3 e 5(11), n = l, 
aiP e 5(12), n<l; t'^^^^t\+2q 
ai(3 e 5(12), n = l, 
aiP e 5(10), 
(ai,/3) = (1,1), n= 1, 
(ai,/3) = (1,1), n^2, 
(ai,/3) = (1,1), n>2, 
61 ^ (0), 62 = (0) 



bi ^ (0) 

fei 7^ (0) 
61 = (0) 
61 = (0) 



^2 = (0) 

b2 = (0) 

^2 ^ (0) 

^2 ^ (0) 



S'fep 6. Define the image of a map R 



1^1, 
I ^ I, n < I, 
I I, n = I, 
l^l, 
1^1. 

Aff{Bi)®Aff{B,) 



j,72 + H(6l®fc2)5' (g, ;j7l--H(fcl 



Aff{B,)(®Aff{Bi) by 



(3.12) R : z'^^bi 

Now the theorem is 

Theorem 3.4. Assume i2.13\) . Then the map \ 3. 12(1 is the combinatorial R matrix 
for U'q{Df^) -crystals. 

Example 3.5. The I — 3 case. 

R : z''i(221)(^z''^(2) ^ z''^+^-'^\l) ® z^^-'^-^\02). 

Because of the bumping algorithm and reverse bumping algorithm, we have 

(221) * (2) = (2 ^ 221) = 021, i7((221) ® (2)) = -2, 

021 = (0 ^ 21) = (0 (2 ^ 1)). 

Example 3.6. The / = 4 case. 

R : (301) ® z^^ (3) ^ z'^^+(-'^) (2) ® z'^i-(-2) (2221) 

Because of the bumping algorithm and reverse bumping algorithm, we have 

3 1 
3 



(301) * (3) 



i?((301) ® (3)) = max (-2, ((301) * (3))i - 4 - 1) = -2, 
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3 1 
3 



1 
2 



2^ ] = {2^[2^ ^ jj =(2-(2-(2-(3^0)))). 

Since (3 — ^ 0) is replaced to (1 — > 2) in / = 4 case, we have 

(2 - (2 -> (2 - (1 - 2)))). 



3 1 
3 



3.3. Combinatorial R for the highest weight elements. We call an element 
& of a ^7^(1^4 crystal a Uq{G2) highest weight element if it satisfies 

(3.13) i,b = fori = 1,2. 

Let 6i (8) ^2 G -B; (8) Bi be a Uq{G2) highest weight element. Then, 

(3.14) hx®b2^ (71,0,0,0,0,0) ® {xi,X2,xz,X3,Q,xi) 

(3.15) (a;2, x^, is, ^i) = (0, 0, 0, 0), 5(63) < 1 for n = 0, 

(3.16) (xa, X3) = (1, 1), (0, 0), 5(62) < 1 if n = 1, 

(3.17) a;3 ^ 2, 5(62) < 1 otherwise. 

For any m G Z>o, we denote (to, 0, 0, 0, 0, 0) £ Bi by (1™). For example, 
(0,0,0,0,0,0) = (0)^(1°), 

(to, 1,1,1,0,0) = (20) if TO = ; = (1"20) otherwise. 

Proposition 3.7. Let a map l : Bi® Bi ^ Bi® Bi he the crystal isomorphism for 
Uq{Df^) crystals. For < n < I set &i «) 62 = (1") ® (/?) e Bi®Bi. Then we have 



(3.18) 



iih ® 62) 



'(1)^ 


5 (l"+il) 




1, 


< 


n<l- 


2, 


i<t>)i 


3(1') 


= 


1, 


n = 


/-l. 




(l)s 


5(1') 


= 


1, 


n — 


/, 




(1)^ 


5 (l"-i20) 


= 


2, 


1 < 


n<l- 


1, 




?)(1'~'2) 


= 


2, 


n ~ 


l, 




(1)5 


5 (l"-iO) 




0, 


1 < 


n<U 




' (1)5 


5 (l"-22) 




3, 


2 < 


n<l. 




(1)5 


5(1) 




1, 


n — 


1, 




(1)5 






1, 


n — 


2, 




(1)5 


5 (l"-2) 


*//3 = 


1, 


3 < 


n<l. 




(0)5 


3 (1") 


*//3 = 


<!>, 


< 


n<l- 


1, 


1(1)5 


5(1'"') 


*//3 = 


<!>, 


n — 


l. 





We will give a proof of this proposition in Appendix A. 



Corollary 3.8. The energy function associated with the crystal isomorphism l\3.18\i 
is given by 



(3.19) 



H{bi ® 62) = 



-2 i/6i®62 = (l")<8(0), 

max(— 2, (bi * 62)1 ^ ^ ^ 1) otherwise. 
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Proof. Thanks to Example 13. 1|, we have 



{bl * 62 



((!")* (/3))i 



if n = 0, /3 = 0, 

1 if n = 0, /3 7^ 0, 
n + 1 if 0,(3=1, 
n if n ^ 0, /3 = 2,0,0, 
n-1 if n ^ 0, /3 = 3,1. 

Thanks to Appendix A, the energy function H{bi 62) is given as foUows. 



Hibi 



'0 


if 1 = 1, bl 


® 62 


= (1)0(1), 


-1 


if 1 = 1, bl 


® b2 


= (1)0(2), (1)0(0), (0)0(1), 





if 1 = 2, bl 


® ^2 


= (11)0(1), 


< -1 


if 1 = 2, bl 


® 1)2 


= (11)0(2), (11)0(0), (1)0(1), 





ifl>3, bl 


® ^2 


= (1')0(1), 


-1 


ifl>S, bl 


® 62 


= (1')0(2), (10 0(0), (l'-i)0(l 


.-2 


otherwise. 







Hence the statement is vahd. 



□ 



3.4. Column insertion and Uq{G2) crystal morphism. For a dominant integral 
weight A, let B{X) be the Uq{G2) crystal associated with the irreducible highest 
weight representation ^(A) [12]. The elements of B{\) can be represented by 
semistandard tableaux for G2 of shape A. 

Let B{\) -B(A') ~ 0^- B{\j)®'^^ be the tensor product decomposition of crys- 
tals, where Aj's are distinct highest weights and mj{> 1) is multiplicity of B(Xj). 
Forgetting the multiplicities we have the canonical morphism from B{X) i3(A') to 
B{\j). The following proposition is due to [17]. 

Proposition 3.9. Let bi * 62 be the tableau obtained from successive column in- 
sertions of letters of the Japanese reading word of 62 into bi. Define a map 
$ : B{X) B{X') — > 0^. B{Xj) by $(61 62) = 61 * 62- Then, the map $ gives the 
unique crystal morphism. 

Thanks to this proposition, we have 

Lemma 3.10. Let l : Bi® Bi ^ Bi® Bi be the crystal isomorphism for U[j{D\ ')- 
crystal. Let 5i 62 G 5/ Bi and 62 b'l € Bi ®) Bi be the Uq{G 2) -highest weight 
elements. For given 61 62 S Bi and 62 b'l € Bi ® Bi, assume that 

(1) bi*b2=b'^*b'l. 

(2) bl 62 (resp. 62 * ^i) i^ connected to bi 62 (resp. 62 b'l) by 1- and 2-arrows. 

(3) i(bi®b2) = b'i®yi. 
Then we have i(5i 62) = b'2 b'l. 

3.5. Proof of Theorem l3.4l First, we show the statement for the energy function. 
Since /i and /2 does not change the shape of the tableaux [12] and H{bi 62) = 
H{bi 62), the algorithm for energy function is reduced to CoroUarv 13.81 Hence 
we have 

f-2 if e 5(10), a27 e S(ll), 

H{bi 62) = H{bi 62) = ; '\ . 

I max(— 2, (61 * 02)1 — t — 1) otherwise. 
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Next, we will show the statement for the crystal isomorphism. By definition, the 
explicit form of b'2 <8) b[ G Bi ^ Bi is given by 



(1)0 


(/?ai ... 


1) 


if ai/3 e B(ll), n < ; - 1, 


(<^)« 


(/3 ai ... 


cei-l) 


if aiP e B(ll), n = l-l, 




g) (/3 ai . . . 


ai-i) 


if ai/3eB{n), n = l, 


(1)® 


(1) 




if (ai,/3) = (1,1), n= 1, 


("2) 


»((/)) 




if (ai,/3) = (1,1), n = 2, 


(a„) 


(g) (a2 0:3 • 


■ ■ Ctn-l) 


if (ai,/3) = (1,1), n> 2, 


(an) 


® (7 a2 as 


■ ■ ■ an-i) 


if ai/3 e 5(10), a27 e S(ll), 


(po)<: 


§> (r;_2 <_ 


-3 • • • r[ qo) 


if ai/3eB(10), a2jGB{12), 






2 ... r[qo p'o) 


if ai/3 G B(12), n<l, 


(Po) C 


3 (^/'-i n'-s 


... r'l qo) 


iiaiP € B{12), n = l, 


(ai) 


»(</-) 




if 61 ^ (0), 62 = ((/)), / = 1, 


(</-)« 


(a\ a2 ■ ■ 


an) 


if 61 ^ (0), 62 = (</-), z 1, n < ; 


(az)c 


g) (ai 0:2 . . 


■ m-i) 


if h^{<t>),b2 = {<f>),l¥=l,n = l 


(<^)8 


(/?) 




ifb, = i4>),b2^i^),l^l, 


(1)0 


(/? 1) 




if fei = ((^), 62 ^ (</-),« 7^ 1, 




(0) 




if 61 = (,^), 62 = (<^). 



where the new letters are defined as follows. 
(i) The aiP e B{10) case. Put 

7 = e(ai/3). 

(ii) The ai/3 € B{10), a2j e B(12) case. Put 

f (r'i,r2, . . . ,r„_2,p„-2,gn-2) = (a„,a„_i, . . . ,a3,Q!2,7), 
\ pi-iqi-irl = r]~'^{riPiqi) for n - 2 > i > 1. 

{Hi) The ai/3 e -B(12), n = Z case. Put 

f (n,r-2, . . . ,ri-i,pi-i,qi-i) = {ai,ai-i, . . . ,a2,ai,/3), 
1 Pi-i9i-ir-- = rj-^iviPiqi) for Z - 1 > z > 1. 

(w) The ai/3 € B(12), n < « case. Put 

(ri,r2, . . . ,r„_i,p„_i,g„_i) = (q;„,q:„_i, . . . ,a2,ai,(3), 
Pi-iqi-irl = V~^{riPiqi) forn - 1 > i > 1, 

PoPo =r'(po). 

Let 5i (g) 62 (resp. 63 ® b[) be the Ug(G2) highest weight elements which arc 
connected to bi 62 (resp. b'2 6'^) by 1- and 2-arrows. By direct computations, 
we have the following equation. 
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'(1") 




if e S(ll), l<n<l, 




(1») 


»(1) 


if (ai,/3) = (1,1), l<n<Z, 




(1") 


»(0) 


if aiP E B{10), a2j E 5(11), 1 < n 




(1") 


»(3) 


if ai(3 E 5(10), a27 e 5(12), 2<n 




' (1") 




if ai/3 G 5(12), 




(1") 




if h 


^ ((/)), 62 - (0), 1 < n < 






)(1) 


if 6i 


= ((/)), 62 (0), 






)(0) 


if bl 


= ((/)), 62 = (0). 








11 Clip t il>l^iij, 1 \ 7T- \ t — Z, 




((h) (X 






if e 5(11), n = ;-l. 






(10 




if e 5(11), n = 






(1) 




if (ai,/3) = (l,i),n=l, 






(0) 




if (ai,/3) = (1,1), n = 2, 




(1) (X 


(1"- 


') 


if (ai,/3) = (l,i),3<n<Z, 






(1"- 


10) 


for ail3 E 5(10), q;27 e 5(11), 


1 < 




(1"- 


^3) 


for ai(3 E 5(10), ^27 E 5(12), 


2 < 




(1"- 


120) 


if e 5(12), 1 < n < ; - 1, 








^2) 


if aiP E 5(12), n = l, 










if 61 ^ (</.), 62 = (<^), ; = 1, 






) (1") 




if 61^(0), 62 -(<^), i<n<; 


-1, 




(1'- 




if 61^(0), 62- ('^),n = ;. 






)(i) 




if 61 = (0), 62 ^ (0) J = 1, 






(11) 




if 61 = (0), 62 ^ (0) J ^ 1, 






)(0) 




if 61 = (0), 62 -(0). 





Thanks to Lemma [3771 and Lemma r3.10[ we have t(6i (8) 62) = 62 ® b[. 
Therefore the proof is finished. 



4. SOLITON CELLULAR AUTOMATON 

4.1. States and time evolution. Let {5/};>i be the family of C/^(5'4'^'')"Crystals. 
Fix a sufficiently large positive integer L. Put u/ = {I, 0, 0, 0, 0, 0) S 5;. We define 
a set of paths by 



(4.1) 



J'l { P = 61 «) 62 «) ■ • ■ 6l e 5f ^; bj = ui if j > 1 } 
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By the iterating combinatorial R: Aff(Bi) (g) Aff(Bi) ~ Aff(Bi) ® A&{Bi), we 
have 

(g) ^°6i O ^°62 O • • • 

H^- z"''h2 (gi ^;-('f^i+'f^2)6(2) ^ ^jOj,^ (g) . . . ^^6^ 

• • • • • • 



i=i 

Here fe*^-^) G is uniquely determined from 6^^^ G i?; and p G CP^ . 

We denote the iterating crystal isomorphism Bi®Bx ^ Ex® Bi by the following 
symbol. 

hi 62 63 



6(0)- 



-6(1)- 



-6(2)- 



-6(3)- 



6(^-1)- 



-6(^) 



6', 



6'. 



Set 6(0) = ui. Assume p = 61 (g) 62 (g) • • • (g) 61, is the element of Ti,. Then, 
b^^^=ui, 6i(g)62 0---(g)6i, gTl, £';(p)<oo. 

Thus we define a map Ti -."Pl — > by 
(4.2) Ti(6i (g) 62 (g) • • • (g) 6l) = 61 (g) 62 (g) • • • (g) 6l. 

We call the dynamical system {^/(p) \ t >0, I > l,p G Jl} the £'4^' -automaton. 
Example 4.1. Using the algorithm given in Theorem 4.4, we have 







1 



1 



1 



111- 



-13- 



■12- 



-21- 



-121- 



-112- 



-111- 



-111- 



1 1 3 1 2 1 1 

Here the frames are omitted. Thus the time evolution is given by 

T3:20O(g)3(gl(g)l(gl(g)l(g--- i-> I0l(g)300(g)l02(g)l(g---. 

Lemma 4.2. For a fixed element ofPL, there exists an integer Iq such that T; = Ti^ 
for any I > Iq. 

Proposition 4.3. (Conservation laws). For any element p G we have 

(4.3) TiTi,{p) = Ti,Ti(p), Ei{Ti,{p)) = Ei(p). 

Proof. Thanks to the Yang-Baxter equation, we can show this by the same argu- 
ment as in Proposition 2.3 of [7]. □ 
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(3) 

4.2. Solitons. In the Dl -automaton a state of the foUowmg form is caUed an 
m-sohton state of length (Zi, ■ ■ ■ , Im) 

(4.4) ...[^i] N ••• M 

Here ...[/]... denotes a local configuration such as 

(4.5) •• -0(1)0(1)® (3)®"^ (2)®"i (1) (1) ® ••• , xi+X2 = l 
sandwiched by sufficiently many (l)'s. 

Remark 4.4. A one-soliton state in the yll^^-automaton is characterized by the 
condition Ei{p) = 1. We also adopt the definition of one-soliton state in D\ - 
automaton from it. 

Lemma 4.5. Let p be a one soliton state of length I, then 

(4.6) The k-th conserved quantity of p is given by Ek{p) = min(fc, /). 

(4.7) Tk{p) is obtained by the rightward shift by min(fc, Z) lattice steps. 

Under the time evolution Tr we identify the local state "... [^] ... " at time t with 

z'^{xi,X2); Xi > 0, X2> 0, Xi+ X2 ^ I, 

(4.8) 

7 — min(r, l)t + (the position of "...[^]..." from the right end). 

Here 7 denotes the phase of "... [i] ..." . From Lemma H75l the phase 7 is invariant 
unless it interacts with other solitons. 

Let Bi (l <E Z>i) be a f7^(A["'^-')-crystal, namely 

(4.9) i3l = {ixi,X2) I Xi > 0, ^2 > 0, Xi+X2= I}. 

Proposition 4.6. Define a map li : Bi U {0} ^ Bf U {0} by 

ii{b) = (3)®"=^ ® (2)®"^! ifb^O ; =0 otherwise. 

Then we have 

iiieiib)) = e2{ii{b)), LiChih)) ^ hMh)) for any b ^ Bi U {Q} . 
Thus each local state . . . [^] . . . is labeled by C/q(Aj^-')-crystal Aff(i3/). 

Remark 4.7. When g„ is a non-exceptional algebra, a one-soliton state is labeled 
by the smaller algebra ?7g(g,i-i). 

Example 4.8. Using the algorithm of combinatorial R given in Theorem 4.4, we 
have r, r, r. 



1 



1 



1 



1 



111- 



-113- 



-133- 



-233- 



-123- 



-112- 



-111- 



1 1 1 3 3 2 1 

Thus the array "332" behaves as a soliton just like in the classical soliton theory. 



We denote a m-soliton state with length {li,l2, ■■■,lm) by 



(4.10) 



'61 



'6™ G Aff(B,J® Aff(B,^ 



•^S{BiJ. 



Assuming li > I2 > 
(4.11) 



> Im, we can expect that the state turns out to be 

...[Ira] •■• [h] [h].... 
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after a sufficiently long time because the longer soliton moves faster under the time 
evolution Tr {r ^ 1). We describe such a scattering of solitons as follows. 

S„, : AS{Bi,) ^ AS{Bi,) ^ ■ ■ ■ ^ AS{BiJ Aff(^;„) ® • ■ • ® Aff(B,J ® Aff(B, J 

The map S2 signifies the two-body scattering of solitons. Now the theorem is 

Theorem 4.9. Assume li > l2- Under the time evolution {r > I2) the map 

S2 is described by the combinatorial R matrix for Uq{A^^^)- crystals. Especially, the 
phase shift is given by 

(4.12) 72 - 72 = 71 - 7i - 2^2 + 3 X H{bi ® 62). 

Remark 4.10. The combinatorial R matrix in Theorem 14.91 has an extra term 2I2 
in the power of z. However, the Yang-Baxter equation holds as it is. 

Corollary 4.11. Scattering of solitons is factorized into two-body scattering. 

The map decomposes into the combinatorial i?'s corresponding to pairwise trans- 
positions of the components. It is independent of the order of the transpositions 
due to the Yang-Baxter equation. For example, we consider the 3-body scattering 
of solitons with length {hjhjs)- Assume li > I2 > Is- Under the time evolution 
Tr {r > li), the rule of 3-body scattering is given by 

(4.13) S3 : z^'bi ® zT^52 ® z'^^'ba ^ z'^'^b'^ ® z'^'^b'^ ® z'^'^b'^ 

The map S3 is factorized into {S2 «) 1)(1 «) S'2)(S'2 1) or (1 ® S'2)(S'2 O 1)(1 » 5*2). 
Since S2 — R satisfies the Yang-Baxter equation, we have 

(^2 <^ 1)(1 ® S2){S2 ® l)(z^^6l ® Z-'^b2 ® Z-f-'b^) 

= (1 ® ^2)(^2 ® 1)(1 ® S2){z^'bi ® Z^^b2 ® Z^^b3). 

4.3. Examples. Let us present examples of soliton scattering. Fix L — 50. For 
a given state p G CP^ compute a series of calculations of ® Bi ~ Bi (g) Br 
by the algorithm in Theorem 13.41 and obtain the time evolution Tr{p). For a 
given two-soliton state z^^{xi,X2) ® z^^{yi,y2) G Aff(i3;J ® Aff(i3i2), put H = 
2I2 — 3min(a;i,y2)- Let 6 be the phase shift. 

Example 4.12. The r — i case. 

t^O: 2222 1111122211111111111111111111111111111111111111 

t^l: 1111 2222 111122211111111111111111111111111111111111 

t = 2: 11111111 2222 11122211111111111111111111111111111111 

t = 3: 11111111111122211122221111111111111111111111111111 

i = 4 : 1111111111111112 2211 1 122 22111111111111111111111111 

S2 = R: z4^(4, 0) ® z38(3^ 0) ^ z38+^(3, 0) ® ^46-^^(4, 0), S = H = 6. 
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Example 4.13. In r = 3 case. 

t = 0: 22211113311111111111111111111111111111111111111111 

t=l: 11122211133111111111111111111111111111111111111111 

t = 2: 11111122211331111111111111111111111111111111111111 

t = 3: 11111111122213311111111111111111111111111111111111 

t = 4: 111111111111 22233 111111111111111111111111111111111 

t = 5: 111111111111111 2203 1111111111111111111111111111111 

t = 6: 11111111111111111123311111111111111111111111111111 

t = 7 : 11111111111111111111130111111111111111111111111111 

t = 8: 11111111111111111111111^31111111111111111111111111 

t = 9: 1111111111111111111111111 11221 11111111111111111111 

i = 10 : 1111111111111111111111111111^21111111111111111111 

t = 11 : 11111111111111111111111111111113321111111111111111 

t = l2: 11111111111111111111111111111111120321111111111111 

i = 13 : 11111111111111111111111111111111111223321111111111 

t = U: 11111111111111111111111111111111111112213321111111 

t = l5: 11111111111111111111111111111111111111122113321111 

S2 = R: z^^3, 0) O ^*i(2, 0) ^ z^'+"{2, 0) ® z^^-"{l, 2), S = H = -2. 

Example 4.14. The r = 3 case. 

t^O: 33211132111311111111111111111111111111111111111111 

t=l: 11133211321131111111111111111111111111111111111111 

t = 2: 11111133213213111111111111111111111111111111111111 

i = 3 : 11111111133232311111111111111111111111111111111111 

t = 4: 11111111111133001111111111111111111111111111111111 

t = 5: 11111111111111131311111111111111111111111111111111 

t = 6: 1111111111111111310O311111111111111111111111111111 

t = 7 : 11111111111111111311133311111111111111111111111111 

t = 8: 11111111111111111131111203311111111111111111111111 

t = 9: 11111111111111111113111112233311111111111111111111 

t = lO: 11111111111111111111311111122133311111111111111111 

i = 11 : 11111111111111111111131111111221133311111111111111 

i = 12 : 11111111111111111111113111111112211133311111111111 

53 : ^^^(1, 2) <g) z^^{l, 1) (g) ^38(0, 1) ^ z^^{0, 1) ® z^\2, 0) z^' (0, 3) 

z^'^il, 2) ® z^\l, 1) ® ^3«(0, 1) ^z^^+^l, 1) ^4^-1(1, 2) ® ^38(0^ 1) 

^^*3(1, 1) ® (1, 0) ® 3) 

^z^'+\0, 1) ® ^^3-2(2, 0) ® ^^^(0, 3) 

^4^(1, 2) ® ^42(1, 1) ® ^38(0^ 1) l®^z47(l, 2) z38+(-l)(i^ 0) ^ z42-(-l)(0, 2) 

i^^39(0, 1) ® ^^^+(-2) (2, 0) ® (0, 3) 
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4.4. Operator Tt,. Let be the crystal given in section [^T^ Put : 
By the iterating crystal isomorphism B^® Bi ^ Bi® Bu we have 



' bL bi u. 



(1) 



>bL 



Here mi, G St, is uniquely determined from p ~ bi 
Define a map Ti, : Bf^ Bf^ by 



i-^ 61 ® • • 

bL e ^L. 



'6l 



(4.14) Ti,(6i (g) 62 (g) ••• ® 6l) = 61 (g) 62 «)•••«) 

Now we identify the two-soliton state 

p = (1)®= ® (3)®"^ (g (2)®"i (g (1)®'^ (g (3)®^^ (g (2)®^i ® (1) ® • • • e J'l, 

(xi + a;2 = /, yi+y2 = k, d > 1) 

with the tensor product of the element of C/q(yl$^^'')-crystal 



(4.15) 



L — c — l — d—k 



(yi,y2) e Aff(BO(gAff(Sfe). 



Lemma 4.15. For the two-soliton state with highest weight z'^^ {I, 0) g) z'^^ (j/i, ^2) G 
^jg^(SO (g ^j9-(Bfe) we /ia«e 



zTi(Z,0)g)zT^(A:,0) */2/2 = 0, 

z^i(/,0)®z^^-3(y^ + l^y2_l) 1/2/2/0. 



Instead of proving the lemma, we give examples below. (See also Appendix B.) 
Example 4.16. Omitting the frames, T\f{z^ {xi,X2)) is described as follows. 



T^(zn4,0)) = 2^4,0). 
1 1 1 



1 



1 1 1 

T^(z^(3,l))-z^(4,0). 



h 



2i- 



1 1 1 

Th(z^(0,4)) = z^(l,3). 
3 3 3 



h 



1 



1 



Lemma 4.17. Assume that I > k. For the two-soliton state with highest weight 
z'f^bi g) z'i^b2 e Ajf{Bi) ® Ajf{Bk) we have 



mR{z^'bi ® z^'b2)) = nmiz^'bi ® z^'b2)). 
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Proof. Put z'^ifoi (8)z'>'=62 = z'^^ 0) ® z'^^ (yi, ya) G Aff(i3/) AfF(Bfc). The y2 = 
case is obvious. Assume ?/2 7^ 0, then we have 

mR{z'''bi(g) z''^b2)) =z'<^+^''-^y^{k,0) ® z'<^-^^+^y^-^{l - y2 + I,y2- 1). 
T^(z''i6i ® z^'b2) ^z^Hl, 0) z^'-^(yi + 1, 2/2 - 1), 

R{Ti^{z^'hl (g) ^'^^62)) =z(T2-3) + (2fe-3(y2-l))(fc^ 0) ® zTl-(2fc-3fe-l))(; -y^ + \, y^- 1). 

Hence we have Tt,i? = i?r^. □ 
Lemma 4.18. mTi{p)) = r,(rn(p)) /or any l>0,pe7L 

Proof. It follows from the Yang-Baxter equation on Aff(i3|,)(g)Aff(i?;)«)Aff(Sp). □ 

Lemma 4.19. For any l^k ^ Z>0; let z^^bi®z''''^b2 and z'^^b'i'^ z'^^b'2 be the highest 
weight elements in Ajf{Bi) (g) Ajf{Bi:). Suppose that they satisfy the conditions 
Tt,(z'>'i5i ® z'>^b2) = mz'^'^b'i (g) z''^b'2) and wt{bi (g) 62) = wt{b[ ® 62). Then we have 
z^^bi ® z^'^b2 = z^'^b\ g) z'^sfea- 

Proo/. Put zTifci (gzT2 62 = zTH^,0) ® z''2(/c- s,s) and z^J 6i g) z'^^ fo^ ^ zTJ(r,0) (g 
z'^^{k' — s',s'). By weight condition to show is in the cases (i) I + k = I' + k', 
(s, s') = (0, 0) or (ii) l + k = l' + k', (s, s') 7^ (0, 0). The case (i) is obvious. In (ii) 
case. By the condition mz'i^bi g) 2^262) = ^{z''''ib[ g) z^^b'^) we have 

z^i(Z,0) g)z^^-3(^ - s + l,s - 1) = z'^i(r,0) g)z'^2-3(/j' _ i^s' _ 1). 

Hence we have (71, 72, ^, k, s) = (71, 72, ^'7 k' , s'). □ 

4.5. Proof of Theorem 14.91 For any two-soliton state z''''^ ^2) ® z'''^ (j/i, ^2) G 
Aff(B/j^) (g Aff(i3/2), Tr (r > 1) commutes with 62, f2- Thus it is enough to check 
the rule for the highest weight elements z'^i 0) (g zT2(t/i, t/2) e Aff(i?; J g) Afr(i3/ J. 
We will show the statement by the induction on 2/2- Because of Proposition [43l we 
have T^T['^_^_i = T^'^^iT^ (r > 12). If a and b are sufficiently large, we can reduce 
the observation of the scattering by Tr in the right hand side to that by Ti^+i in 
the left hand side. Hence the statement is reduced to T/^+i. 
The 2/2 = case. By direct calculation we obtain 

(4.16) 5*2 : z'^i {I, 0) (g z'^^k, 0) ^ z'^^+^'^ik, 0) (g z'^^-^''{l, 0). 
The 2/2 > case. For the two-soliton state 

p =(i)«5(i-'i-7i) ^ (2)»'i (i)®(7i-72~/2) ^ (3)»!y2 ^ (2)«^«i ® (1)^^^ e Tl 
=zTnZi,0) ® z^- (2/1, 2/2) e Aff(SjJ ® Aff(B,J, 
the statement to show is the following: 

(4.17) T^{p)^R{p) forr>/2,i>l. 

Thanks to Lemma HTTl and Lemma HTTSl we have (i?T[,)(p) = {T^R){p) and 
{T^Tr){p) — (T,*Tt,)(p) respectively. By the induction hypothesis we have (i?T[,)(p) = 
{T^^)Ip). Hence we have T^{T^{p)) = T^{R{p)). Since both r,*(]5) and i?(p) have 
the same weight as p, we can apply Lemma 14.191 Thus we have (|4.17p . 
Therefore the proof is finished. 
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Appendix A. Proof of Proposition 13.71 
Because of the weight, we have 

(A.l) L : {1, 0, 0, 0, 0, 0) ® (1, 0, 0, 0, 0, 0) ^ (1, 0, 0, 0, 0, 0) ® (/, 0, 0, 0, 0, 0) 



under the crystal isomorphism l : Bi ® Bi ^ Bi ® Bi. 

Set u„ = (n, 0,0, 0,0, 0) G Bi for < n < L For any 6i (8) &2 G ^ Bi there 

exist gii, - ■ ■ ,gi^ G {e^, fi \i = 0, 1, 2} such that &i (8 62 = 9ii ■ ■ ■ ^ ui). Set 

&2 (8) &i = 6(&i (8) 62)- Then &2 8) /^i = 5ii ■ • -Si^Cui 8) ui). Denote it by 61 8) 62 — 

h ® • • - Si™- 

By using the foUowing data: 



(A.2) 



/P(0, 0,0, 0,0,0 



(A.3) /^'(0,0,1,1,0,/-1) = 



(A.4) /^'(0, 0,0, 2,0, /-I) 



(A.5) /^'(0,1,0,0,0,/-1) = 



(A.6) /^'(l, 0,0, 0,0, /-I) 



(0,0,0,0,0,/-p) if0<p<Z, 

(p- 00,0,0,0,0) if? + l<p<2L 

(0,0,1,1,0,/- 1 -p) ifO<p</-l, 

+ 1,0,1,1,0,0) if/<p<2/-2. 

(0,0,0,2,0,/- 1 -p) ifO<p</-l, 

(p-/, 1,0, 0,0,0) if/<p<2/-l. 

(0,1,0,0,0,/- 1 -p) if0<p</-2, 

(p-/ + 1,1,1,1,0,0) if /-I <p<2/-3. 

(1,0,0,0,0,/- 1 -p) if0<p</-2, 

(p-/ + 3,0,0,0,0,1) if / - 1 <p < 2/ -4. 



we obtain the crystal isomorphism 61 8) 62 — &2 8> 61; • • • gi„ 
The / = 1 case. 



(1) ® (</)) ~ (1) ® (I); e~o, 
(1)® (1)2.(1)0(1); e~2. 



((/.)8)(0) ^ (0)® (0); eo-B^gi], 
(l)®(0):i^(l)®(0); fofihfiE^el 



(4>) ® (1) ~ (</)) ® (1); i;ig3, (1) ® (2) ^ (1) ® (2); fofifofif2.fiE^el 



The / = 2 case. 



(11) ® (</.) (1) ® (1); go, 
(11)^(1)^(1)0(0); g2, 
(1)®(T)=. (1)0(1); g3. 



(</))®(l)~ (1)0(11); i?ig4, 
(1) 8. (2) ~ (1) ® (20); fofiE'it 
(1)®(0)~(1)®(0); fofihfiE^it 



(1) ® (1) (</.) ® (11); E^el (11) ® (0) ^ (1) ® (10); f^fME^'el 
(1) ® (,/.) ~ (0) ® (1); e~oii;2e5, (n) ^ (3) ~ (1) ® (2); fJJofihfiE'Sl 
{4>) ® {4>) ~ (0) ® (0); e~2i?2g5, (11) ^ (2) ^ (1) ® (12); fJofif2fiE^el 
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The I > 3 case. 

(1'-") ® (1) ~ (1) ® for < n < Z - 2, 

(1)0(1)^(1)0(1); e~[,+\ 
(0)0(1)~(1)0(1I); i?ie~[,+2, 

(1") (2) ~ (1) ® (l"-i20); fSfiE^e^f^^ for 1 < ti < ^ - 1, 

(1") (0) ~ (1) (l"-iO); /r'ei/o/f /2/ii?^eJ,+2 for 1 < n < /, 

(1") (3) ~ (1) ® (l"-22); jSjfhhE^~e'+^ for 2 < n < / - 1, 

(10 (2) ^ (1) ® (l'-i2); fl+'flhhE%+' 

(1") (1) ~ (1) ® (l"+il); for 1 < n < ; - 2, 

(i'-i)®(i)^(</')®(iO; i^'^r', 

(1'-") ® (</>) ~ (0) (1'-"); eg£;'eo'+^ for 1 < n < L 
Here i;" = e^e^^ef" S^e^ . 



Appendix B. The table of B^® Bi ~ Bi® B^ 

The foUowing data is all images of the crystal isomorphism B\^ (g) Bi ^ Bi (g) B^ 
that is needed for Lemma [4. 151 

*(12)(g)(l) ~ (l)(g)*(12) *(12)®(2)~(2)®*(12) *(12)(g)(3) ~ (l)(g)(l)2 

(1)2 (8) (1) ~ (1) ® (l)i (l)i®(l)~(l)®*(23) *(23)®(1)~(2)®*(13) 

*(13)(8)(1) ~ (1)(»*(13) (1)2 0(2) ~ (l)®(2)i, (2)i(8(l)~(l)(g)*(20) 

*(20)(8)(1) ~ (2)®*(23) (2)i(8)(2)~(l)(8)*(23) *(23) (g) (2) ~ (2) *(23) 

*(23)(g(l)~(2)®*(20) (1)2® (3) ~ (l)®(3)i (3)i(g(l)~(l)®*(30) 

*(30)(8)(1) ~ (3)(»*(23) (3)i®(2)~(l)(8)*(33) *(33) ® (1) ~ (3) *(20) 

*(33)(g)(2) ~ (3)(g)*(23) (3)1 O (3) ~ (1) (8)*(32) *(32) (8) (1) ~ (3) (g) *(30) 

*(32) (g) (3) ~ (3) ® *(32) *(32) ® (2) ~ (3) ® *(33). 

where *(«/?) denotes 
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